Classical Solutions for a nonlinear Fokker-Planck 
equation arising in Computational Neuroscience 



Jose A. Carrillo 1 , Maria d. M. Gonzalez 2 , Maria P. Gualdani 3 and Maria E. Schonbek 4 

1 Institucio Catalana de Recerca i Estudis Avangats and Departament de Matematiques 
Universitat Autonoma de Barcelona, E-08193 Bellaterra, Spain 

2 ETSEIB - Departament de Matematica Aplicada I 
Universitat Politecnica de Catalunya, E-08028 Barcelona, Spain 

3 Department of Mathematics 
UT Austin, Austin, TX 78712, USA 

4 Department of Mathematics 

UC Santa Cruz, Santa Cruz, CA 95064, USA 



Abstract 

In this paper we analyze the global existence of classical solutions to the initial boundary- 
value problem for a nonlinear parabolic equation describing the collective behavior of an ensem- 
ble of neurons. These equations were obtained as a diffusive approximation of the mean-field 
limit of a stochastic differential equation system. The resulting Fokker-Planck equation presents 
a nonlinearity in the coefficients depending on the probability flux through the boundary. We 
show by an appropriate change of variables that this parabolic equation with nonlinear bound- 
ary conditions can be transformed into a non standard Stefan-like free boundary problem with 
a source term given by a delta function. We prove that there are global classical solutions for 
inhibitory neural networks, while for excitatory networks we give local well-posedness of classical 
solutions together with a blow up criterium. Finally, we will also study the spectrum for the 
linear problem corresponding to uncoupled networks and its relation to Poincare inequalities for 
studying their asymptotic behavior. 

1 Introduction 

The basic models for the collective behavior of large ensemble of interacting neurons are based 
on systems of stochastic differential equations. Each subsystem describes an individual neuron 
in the network as an electric circuit model with a choice of parameters such as the membrane 
potential v, the conductances, the proportion of open ion channels and their type. The individual 
description of each neuron includes an stochastic current due to the spike events produced by 
other neurons at the network received through the presynaptic connections. We refer to the 
classical references [HI [TUJ [23] and the nice brief introduction [TSJ for a wider overview of this 
area and further references. As a result of the coupling network, the collective behavior of the 
stochastic differential system can lead to complicated dynamics: several stationary states with 
different stability properties and bifurcations, synchronization, and so on, see [TJ 1171 118] for 
instance. 

To understand this behavior the evolution in time of the potential through the cell membrane 
v(t) has been modeled by several authors [TJ [H [201 El H21 (HI- The neurons relax towards their 
resting potential vl (leak potential) in the absence of any interaction. All the interactions of 
the neuron within the network are modeled by an incoming presynaptic current I(t) given by 
an stochastic process to be specified below. Therefore, the evolution of the membrane potential 
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is assumed to follow the equation 



dv 

C m — = -g L (v~v L )+I(t), (1.1) 

where C m is the capacitance of the membrane and, ql the leak conductance. If the voltage 
achieves the so-called threshold voltage (maximum voltage), the neuron voltage is instanta- 
neously reset to a fixed voltage vr. At each reset time the process produces a spike, which 
builds up the incoming presynaptic current I(t) and is added to the mean firing rate produced 
by the network N(t) defined as the average number of spikes per unit time produced in the 
network. 

Most of the microscopic models for neuron dynamics assume that the spike appearance 
times in the network follow an independent discrete Poisson process with constant probability 
of emitting a spike per unit time v. We will assume that there are two types of neurons: 
inhibitory and excitatory, and each produce a spike of strength Je and Jj respectively at 
their spike times. The total presynaptic current I(t) in ([1.1)1 . coming from the spikes within 
the network, is computed as the difference of the total spike strengths received through the 
synapsis by a neuron at the network composed by Ce excitatory and Ci inhibitory neurons. 
This stochastic process I{t) has mean given by fie = Bv with B = CeJe — CiJi and, and 
variance = (CeJe + CiJj)v. We will say that the network is excitatory if B > (inhibitory 
respectively if B < ). Dealing with these discrete Poisson processes can be difficult and thus, 
an approximation was proposed in the literature. This approximation consists in substituting 
the stochastic process I(t) by a standard drift-diffusion process with the same mean and variance 

J(t) dt ps lie dt + a c dW t 

where Wt is the standard Brownian motion. We refer for more details of this approximation to 
[TJE1I20, 6, 22, [19l [15] . The approximation to the original Leaky Integratc&Firc neuron model 
(|l.lj) is then given by 

dv=(-v + v L +/j,c)dt + a c dWt (1.2) 

where we choose the units such that C m = = 1, for v < Vth with the jump process: v(t^) = vr 
whenever at t a the voltage achieves the threshold value v(t~) = v t h] with vl < vr < v t h- The 
last ingredient of the model is given by the probability of firing per unit time of the Poissonian 
spike train v, i.e., the so-called total firing rate. The firing rate depends on the activity of the 
network and some external stimuli, it is given by v = v ext + N(t) where N(t) is the mean firing 
rate produced by the network and v ext > is the external firing rate. The value of N(i) is then 
computed as the flux of neurons across the threshold or firing voltage Vth- 

Studying the stochastic problem (jl.2|l with the jump process specified above can be written 
in terms of a partial differential equation for the evolution of the probability density p(v, t) > 
of finding neurons at a voltage v £ (— oo, vth] at a time t > 0. This PDE has the structure of a 
backward Kolmogorov or Fokkcr-Planck equation with sources and given by 

^(v,t) = — [(v-v L - f i c )p(v,t)}+^-^(v,t)+N(t)6 v = VR , v<v th . (1.3) 

A Delta Dirac source term in the right-hand side appears due to the firing at time t > for 
neurons whose voltage is immediately reset to vr. Imposing the condition that no neuron should 
have the firing voltage due to their instantaneous discharge, we complement (|1.3j) with Dirichlet 
and initial boundary conditions 

p(v th ,t)=0, p(-oo,t)=0, p(v,Q)= Pl (v)>0. (1.4) 

The mean firing rate N(t) is implicitly given by 

N(t):=-^f^(v th ,t)>0, (1.5) 
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that is the flux of probability of neuron's voltage, is the least at v t h- It is easy to check that 
this definition implies, at least formally, that the evolution of (|1.3j) is a probability density for 
all times, that is 

p(v,t)dv= / pi(v)dv = l 



oo 



for all t > 0. Let us note that in most of the computational neuroscience literature [TJ [15] . 
equation (|1.3|) is specified on the intervals (—00, vr) or (vr, v t h) with no source term but rather 
a boundary condition relating the values of the fluxes from the right and the left at v — vr. The 
formulation presented here is equivalent and more suitable for mathematical treatment. Other 
more complicated microscopic models including the conductance and leading to kinetic-like 
Fokkcr-Planck equations have been studied recently, see [4] and the references therein. 
Finally, the nonlinear Fokker-Planck equation can be rewritten as 

dv (T 2 d 2 v d 

^ = Y ^ + g- v [(v-fi)P} + N(t)S^ R ,v< Vth 

where a 2 = 2a 2 , + aiN(t), with ao > 0, a\ > and jl = Bv ext + BN(t). We will focus only on 
the simplest case in which the nonlinearity in the diffusion coefficient is neglected by assuming 
a\ = 0. Without loss of generality, we can choose a new voltage variable v < and an scaled 
density p defined by 

p(t,v) = f3p(t, f3v + vth) 
where /? = ao- Then our main equation, after dropping the tildes, reads 



% = Pt + £-[(V-l*)p}+X®tv=v H , v<0, (1.6) 
at ov ov 



where the drift term, source of the nonlinearity, is given by 

fi = b + bN(t) withJV(t) = -^(0,t)>0 (1.7) 

ov 

with bo = (Bv ex t — vth) I ao and b = B/clq. Let us remark that the sign of bo determines if 
the neurons due only to external stimuli may produce a spike or not, therefore it controls the 
strength of the external stimuli. 

In a recent work [3], it was shown that the problem (|1.6[) - (|1.4[) - (|1.7|) can lead to finite-time 
blow up of solutions for excitatory networks b > and for initial data concentrated close enough 
to the threshold voltage. Here, we give a characterization of the maximal time of existence of 
the classical solution, if finite, and thus, of the blow up time. We show that if the maxima 
existence time is finite, it coincides with the time in which the firing rate N(t) diverges. This 
divergence in finite time of the firing rate has no clear biological significance. It could mean that 
some sort of synchronization of the whole network happens, see [3] for a deeper discussion. This 
is an scenario that does not show up in the typical reported applications [HE]. In the rest of 
this work, we concentrate in studying the existence of classical solutions to the initial boundary 
value problem (|1.6|) - (|1.4[) - (|1.7p . We show that solutions exists globally in time for inhibitory 
networks b < and, we give a characterization of the blow up time for the case when b > 0. 
Although the precise notion of classical solution will be discussed in the next section, the main 
theorem of this work can be summarized as follows. 

Theorem 1.1. Let pi{x) be a non-negative C 1 ((— 00, Vth]) function such that pi(vth) = 0. 
Suppose that pi, (pi) x decay at —00, then there exists a unique classical solution to the problem 
([j ^ -([P | -([T77 ) on the time interval [0,T*) with T* = 00 for b < and, for b > 0, T* > can 
be characterized by 

T* = sup{i > : N(t) < 00} . 

Furthermore, for b > there exist classical solutions blowing up in finite time, and thus, with 
diverging mean firing rate in finite time. 
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Let us remark that the last statement is merely obtained by combining the result in [3] with 
our classical solutions existence result and the characterization of the maximal time of existence. 

The main strategy of the proof as is shown in section 2, is given by an equivalence. This 
equivalence, through an explicit time-space change of variables, transforms our problem into a 
Stefan-like free boundary problem with Delta Dirac source terms, resembling price-formation 
models studied in jllj . In section 3, we will use ideas and arguments in Stefan problems [8] to 
show local existence of a solution. Next, in section 4 we will prove global existence of classical 
solutions for inhibitory networks (b < 0) and give a characterization of the blow up time for 
excitatory networks (b > 0). The difference between the cases b < and b > corresponds to 
the well studied Stefan-problem in the normal and in the undercooled cases, see [TB] for classical 
references in the Stefan problem. The final section is devoted to study the spectrum of the 
linear version of (|1.6p (6 = 0) that has some interesting features and properties connected to 
classical Fokkcr-Planck equations. 



2 Relation to the Stefan problem 



The main aim of this section is to rewrite equation (|1.6[) as a free boundary Stefan problem with 
a nonstandard right hand side. For this we recall a well known change of variables, [S], that 
transforms Fokker-Planck type equations into a non-homogeneous heat equation. This change 
of variables is given by 

1 



y = e v, t 



that yields 
or equivalently 



(e«-l), 



p{v,t)=e t w[e t v,-(e' t -l)\, 



w(y,r) = (2t + 1)-^ 2 p 



= ,-log(2r + i; 



,\/2T + T 2 

In the sequel, to simplify the notation, we use a{r) ~ (2r + l) -1 / 2 = e _t . A straightforward 
computation gives for w: 



w r = w yy - fi(T)a(r)w y + M(t)5 ' *r 



(2.1) 



where M(r) = a 2 (T)N(t) 
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The additional change of variables: 



u(x, t) = w(y, t) wherex^y— / n(s)a(s) ds = y — bo (Vl + 2f — l) — b / M(s)a 1 (s)d 



removes the term with w y in ()2.1[) . For boundary conditions at initial time, denote s/ = Vth{= 0). 
We have the following equivalent equation 

Lemma 2.1. System (|1.6I) - (|1.4[) - (|1.5[) is equivalent to the following problem 

u t = u xx + M(t)6 x=si ( t ), x <s(t),t>0, 

s(t) = sj - b + 2t - 1) - b [ M(s)a- 1 (s)ds, t > 0, 



si(t) = s(t) 



du 

M ® = " fa 



x— s(t) 

u(-oo,t) = 0, u(s(t),t) = 0, 
u(x,0) = m(x), 



t > 0, 

t > 0, 

t > 0, 
x < sj. 



(2.2) 
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Proof. The proof is straightforward by the changes of variables specified above and, as such is 
omitted. □ 

We now give a definition the concept of classical solution. In what follows we work with 
the Stefan- like free boundary problem (|2 . 2[) . It is immediate to translate this to a concept of 
classical solution to the original problem (|1.6[) - (|1.4[l - (|1.5[l by substituting u by p, x by v, M{t) 
by N(t), si(t) by vr, and s(t) by v t h- 

Definition 2.2. Let uj(x) be a non-negative C 1 ((— oo,s/]) function such that uj(sj) = 0. Sup- 
pose that ui,(uj) x decay at — oo. We say that (u(x,t), s(t)) is a solution of (|2.2[) with initial 
data uj(x) on the time interval J = [0, T) or J — [0, T], for a given < T < oo, if: 

1. M(t) is a continuous function for all t G J, 

2. u is continuous in the region {(x,t) : — oo < x < s(t),t € J}, 

3. u xx and Ut are continuous in the region {(x,t) : — oo < x < Si(t),t G >A{0}} U {(x,t) : 

si(t) < x < s(t),te J\{0}}, 

4- u x (si(t)~ ,t), u x (si(t) + , t), u x (s(t)~ ,t) are well defined, 

5. u x decays at — oo, 

6. Equations (|2.2p are satisfied. 



The next lemma presents some of the a priori properties of the solution to 
Lemma 2.3. Let u(x,t) be a solution to (|2.2[) in the sense of Definition \2.'2\ Then 
i) The mass is conserved, 

s(t) pSl 

u(x,t)dx= / uj(x)dx, 



for all t > 0. 

ii) The flux across the free boundary si is exactly the strength of the source term: 

M(t) := -u x (s(t),t) = u x ( Sl (t)-,t) - u x ( Sl (t)+,t). 

Hi) If bo < and b < (resp. bo > and b > 0), the free boundary s(t) is a monotone 
increasing (resp. decreasing) function of time. 

Proof, i) Mass conservation, follows by integration of the equation and straightforward integra- 
tion by parts. 

ii) To establish the jump across the free boundary, i.e. part ii), integrate the first equation in 
(|2.2[) over the interval (— oo, s±(t)) , yielding 

si(t) /-si(t) 

u t dx — / u xx dx = 0. 



oo 



Hence, 

d r iW 

— / u(x,t)dx = u x (s 1 (t)-,t) + s 1 (t)u{s 1 {t),t). (2.3) 

Ot J-vo 

Similarly, an integration of the first equation in (|2.2[) in the interval (s\(t), s(t)) gives 
d f s{t) 

— / u(x,t)dx + s 1 (t)«(si(t),*) - s(t)u(s(t),t) = u x (s(t),t) - ■u ;c (s 1 (<) + ,t). 
<Jt J Sl ( t ) 

If we substitute u(s(t),t) = in the previous line it follows 



d r s{t) 

— I udx + s 1 (t)u(s 1 (t),t)=u x (s(t),t)-u x (s 1 (t) + ,t). (2.4) 

81 (t) 



dt 
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d f s{t) 

= — / u(x,t)dx = u x (si(t) ,t) +u x (s(t),t) - u x (si(t) + ,t). 



Adding (|2.3[) to (|2.4p and recalling that the mass is preserved we get 

at _ 

It follows that 

u x (s(t),t) = u x (si(t) + ,t) - u x (si(t)~,t), 

as desired. 

Let us prove the last part Hi). The free boundary is a monotone increasing function of time 
since b < 0, b < 0, a > 0, and 



s(t) = S j - b (VI + 2t - l) - b / A/^oT 1 ^)^, i>0, 

Jo 

while the fact that M(t) is strictly positive follows by the classical Hopf 's lemma. □ 



3 Local existence and uniqueness 

In this section we prove local existence of solution. Our method is inspired by the theory 
developed by Friedman in [7J 18] for the Stefan problem. We first derive an integral formulation 
for the problem. A derivative with respect to x yields an integral equation for the flux M, 
where a fixed point argument can be used to obtain short time existence. Once M is known the 
equation for u decouples and it is solved as a linear equation. 

Theorem 3.1. Let uj(x) be a non-negative C 1 ((— oo,s/]) function such that ui(si) = 0, and 
Mi, (ui) x decay at — oo. Then there exists a time T > 0, and a unique solution (u(x,t), s(t)) of 
the equation (|2.2|) in the sense of Definition [2~2l for t £ [0,T] with initial data uj. Moreover, 
the existence time T is an inversely proportional function of 

sup K(x)| . 

— OC<X<Sj 

The proof of Theorem 13.11 will be divided in several steps. The first step deals with an 
integral formulation of the solution, which is used to show the existence of M. 

3.1 The integral formulation 

Let G be the the Green's function for the heat equation on the real line: 

1 f |x-^| 2 

G(x, t, £, t) = -— —-J- exp 



[4?r(t - r)]V2 ^\ 4(t-, 
To obtain an integral formulation of the solution u of (|2.2p , recall the following Green's identity 



To recover u we first integrate the identity (|3.ip in the two regions 

-oo < £ < si(r), < t < i, and Si(t) < £ < s(r), < r < i, 

and then add up the results from the integration. We split the resulting expression into the 
following four terms; the only problematic one is the one containing u^: 

t r»(r) q / QQ\ ft f s(r) g 
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Each term will be analyzed separately. Note that u and G have enough decay as |£| — > oo to 
justify the following computations due to Definition 12.21 Since G(x,t, — co,r) = it holds 



£=«i(t) 



dr 



^ — — OO 



G(x,t,s 1 (r),r) ^ 



dr. 



i(t)" 



Next, we obtain 
v 9u 



7J 



6* 



9* 



- G — 

e=.(r) 9 ? 



€=«i(t)+. 



<^ =s(t) M(r) + G 



du 



(3.2) 



dr. (3.3) 



S= Si (t)+ 



Here we have used that ^| 



7/7 



£= S (t) 

dG 



■M(t). For the third integral we have 
<9G N 



«=«(t) 



dr 



£= — OO 



{{ KS(T) ' T) |f 



- u(-oo,r) — 



(3.4) 



dr = 0, 



^ — — OO I 

taking into account that u(s(r),r) = u(— oo,r) = 0. Finally, using that u(s(r),r) = 0, we have 



iv=ri 



s(r) 



Gud^dr 



lo dT J- 

Recall that G(x,t,£,t) = 5 X =£, thus the last identity yields 



r a(t) r s(0) 

/ Gu\ T=t dt- Gu\ T=0 

J — oo J — OO 



g(f) ,8(0) 

/ G(x,t,^0) Ul (Odt 



(3.5) 



Combining ()3.2|) . ((3T3J) , (|3.4j) . (|3.5[) . and part of Lemma l2~3l yields that the solution u reads 
as 



u(x, t) 



s(0) f* du 

G(x,t,Z,0)MtW+ / G(x,t,«i(r)) — 

-oo JO 

M{T)G{x,t,s{T),T)dT - [ G(x,t,8l(T)) 

Jo 



C=«i(t)- 

du 



(0) f t ,-t 

G(:r,i,£,0)u/(Od£- / M(t)G(M, s(r), r)dr + / M(t)G(x, t, Si(r), r)dr 
» Jo Jo 

= : h-I 2 +h. (3.6) 
The term Ji represents the solution to the heat equation with data u/. Indeed, 

>(0) 

G(x,t,e,0) U/ (0^ 

OO 

is the solution to the homogeneous heat equation with initial data 

„ (n _ / MO e < *(o), 

Uo(c) - \ o c > *(o). 

All the calculations up to here are formal assuming that u is a solution of the equation (|2.2j) 
as in Definition ^. 21 We now derive an equation for M which will be solved for short time using a 
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fixed point argument. The first step is to obtain the space derivatives of the terms i = 1,2,3 
and evaluate them at x = s(t)~: 



dh 

dx 



S (ty 



,■((» 



G x (x,t,Z,0)ui(Od£ 



.(0) 



g(m,£,oK(£K. 



To get the derivative of h, we use [7J Lemma 1, pag 217]: this lemma states that for any 
continuous function p, 



hm — . 

x->s(t)~ OX J 



dG 

p(T)G(x,t,s(r),T)dr = ^p(t) + / p(r) — t, s(r), r) dr 



(3.7) 



As a consequence, 



dh 

dx 



x—s(t)~ 



M{t) + M( T )G x (s(t),t,s(T),T)dT. 



For the derivative of h note that problems can only occur if t = r and s(t) = Si(r), but this is 
not possible by the definition of s± . Thus, 



dx 



x—s(t) ~ 



G x (a(t),t; 3l (T),T)M(T)dT. 



Substituting the estimates on 1\ , h and ^3 into (|3.6[) we get 



S (0) 



-M(t) 



After reordering yields, 

ns(0) 



G(8(t)M,0M{t)d£--M(t) 



M{T)G x {s{t),t,s{r),T)dr+ / M( T )G x (s(t), t; si(r), T )dr. 



M(t) 



G( s (t),t,^0)u'j(Od^ 



+ 2/ M(T)G x { S {t),t, S (T),T)dT-2 M{T)G x { S (t),t, Sl {T),T)dT. 



(3.8) 



3.2 Local existence and uniqueness for M 

Theorem 3.2. Letuj(x) be a non-negative C 1 ^— 00, sj]) function such thatuj(sj) = 0. Suppose 
Ui, (uj) x decay to zero as x — > — 00. Then there exists a time T > such that M(t) defined by 
the integral formulation (|3.8[) exists fort S [0,T] and is unique in C([0,T]). The existence time 
T satisfies 

T < I sup |w.j(ai)| 

\— OC<X<S] 

Proof. The local in time existence of M(t) is obtained via a fixed point argument. For this, we 
modify the classical argument for the Stefan problem to account for the additional source term 
given by M(t) S x=ai uy For er, m > 0, consider the norm 

||M|| := sup \M(t)\ 

0<t<a 



in the space 



Grj 



{M eC{[0,a}) : ||M|| < to}. 
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Set 

, S (0) 



T(M)(t) := - 2 / G(a(t),t,t,OWAt)dS 

J —oo 

(3.9) 



t pt 



+ 2 M(T)G x (s(t),t,s(r),T)dT-2 M( T )G x (s(t),t, si(r),r)dr 
./o Jo 

:= Ji + J 2 + Ja- 

In order to apply fixed point arguments, it is necessary to show that for sufficiently small cr we 
have: T : G a , m — > C a ^ m and that T is a contraction. Define 

m:=l + 2 sup (3.10) 

— oo<a;<s(0) 

5iep L- We show that for cr sufficiently small T : C a _ m — > C G ^ n . For simplicity, we focus on 
the proof in the case b < 0. At the end we make the necessary corrections for b > 0. Choose a 
sufficiently small so that 

i. a _1 (t) < 2, V< < a , 
it. "'d f 'ol+ 2TO l b l) t7 i/2 < U2 

Hi. \vr\ - \b \cr > 0, 
2m f°° 

iv. -= z- 1 exp{-z 2 }dz< 1/2. 

s/8ct 

We obtain first an auxiliary estimate. Since a has been chosen so small that condition i. holds 
and a _1 \/l + 2t is a 1-Lipschitz function for t > 0, if M € C CTTO then 

[»(*) — «(t)| < |6o||*-r| + |6| / M(s)a" 1 (s) ds < (|& | + 2|&|m)|t-r|, (3.11) 



i.e., s(£) is a Lipschitz continuous function of time. 

To estimate the image of the operator T(M) as defined in (|3.9|) for M £ C<j yTn we find 
separately a bound for each term J±, J 2 , J3. First, for Ji, note that 

s(0) 

G(s(t),t,Z,0)d£ < 1. 



Then, it is straightforward to check 



r ) mo) 

|Ji|<2<^ sup \uj(x)\\ G(x,t,£,0)d£<2 sup |u/(a;)|. 



We bound J 2 as 

rt 



\M < 2m / \G x (s(t),t,s(T),T)\dT. 
Jo 

Substituting 



2v^(i-r) 3 / 2 4(t-r) 
in (|3.1ip and taking into account the choice of a given by ii., it follows that: 

m [t \ s (t)-s(T)\ „ f Kt)- a (r)l a ^ 
|J2| -V^X (*-r)3/ 2 eXP |" 4(f-r) / * 



m(|6o|+2m|6|) f* 1 , 2m(|6 | + 2m|6|) 1/2 2m(|6 | + 2m|6|) 1/2 1 

<■ 7= / -, . , /n clt = 7= 1 1 < == a 1 < — 

Jo (t - t) 1 ' 2 v/i^r ~ 2 
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Before we consider J3, we need the following auxiliary estimates. The inequality ye v < e 2 
implies that 

|G *«' T)lfi v5(b) exp {-^}' (3 - 12) 

The definitions of s(t) and si(t) = s(t) + v^a" 1 (r), using that b < 0, and by the condition Hi, 
yield 

|s(t)-ai(T)| > \v R \-\b \a>0. (3.13) 

If we integrate (|3.12[) we get 

|Ga;(s(t),i,si(T),r)|dT < -7= / cxp <^ — r S dr 



v 7 ^ Jo t ~ t r i 8(i - r) 

<^fj- BP Lw-iyu (3.i4) 

~~ V4n Jo t — r \ 8(t - r) J l ° J 

1 f 00 1 2 1 r 1 2 

y^7T / 1 " I — I b o I CT z \pK jknjH^ojf; z 

where we used the change of variables z = \ VR \~\ ha \ a gy the last estimate and by condition iv 

* 2m f 00 1 .2 , 1 



|J 3 |<2m/ \G x (s(t),t, Sl (T),T)\dT<-= - e -»dz<-. (3.15) 

Jo v 71 " J l"«l~ |b ° l ' y 2 2 

The estimates for J^, i = 1,2, 3 establish that T(M) £ C a ^ m since 

1 1 T(M) || < Ji + J 2 + J 3 < m, VM € C a , m , 

by the choice of m in (| 3 . 1 [) . 

It remains to consider the case 6 > 0. It is clear that the only modification needed is the 
estimtc (|3.13p . For this use 

\s(t)- Sl (r)\ = Is^-sCr)-^" 1 ^)! > IMctV) - \s(t) - s(r)\\ > \ VR \-(\b \+m)a, (3.16) 

which may be estimated from below by a positive constant for some a small enough. Then, we 
have the same result as in the case b < assuming analogous conditions to i., ii., Hi., iv. above. 
The main difference between the cases bleqO and b > is that in the case b > in (|3 . 13[) the 
difference between the free boundary s(t) and the source Si(r) for < r < t now depends on 
the bound of the initial data (|3.10|) . 

Step 2.- The mapping T : C a ^ m —> C a ^ m defined in (|3.9[) is a contraction for a small enough. 
In the sequel the constant C is arbitrary and may change from line to line. Let M, M £ C CTim , 
and 



(t) = sj -b (yTT2t- l) -b / M(T)aT 1 {T)dT, . (3.17) 



s{t) =si-b (Vl + 2i - 1) - b I M(r)a _1 (r) dr . 
The following auxiliary estimate holds: 

|s(t) - s(t)\ < \b\ [ \M{r)-M{T)\a- l {T)dT<\b\\\M-M\\ I \J2t + Idr 



^-\\M-M\\ [(2t+l) 3 / 2 -l 
3 



(3.18) 
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It is straightforward from (|3.17p that 

\s(t) - S(t)\ < 2 \b\ || M -M\\, < t < a < 1. 
By condition i. on a and (|3 . 1 1 1) it follows that 

max{| S (t) - s(t)\, \s(t) - s(r)|} < (|6 | + 2m |b|)|t - r| < (|6 | + 2 |fc|)ro|t - r| 
To show that T is a contraction we proceed as follows. 

r s(0) 



(3.19) 
(3.20) 



\T(M)-T(M)\ <2 



|« / / (OI|G(a(t),«,f,0)-G(S(t) > t,f > 0)|de 



M(r)G,(s(i), t, s(r), r) - M(r)G a (S(*), t, S(t), r) dr 



M(t)G x (sit), t, Sl (r), r) - M(r)Gf a (S(t), t, Si (r), r) dr 
= : Ai + A2 +A3. 

Without loss of generality assume that s(t) > s(t). The mean value theorem applied to the 
kernel G(x,t,£,0) gives for some s G [s(i),s(i)] 



Recall that 



\G(s(t),t, £, 0) - G(5(t), t, e, 0)| < \G x (s, t, £, 0)| • |«(t) - a(t)|. 



(3.21) 



iG s (g,t,e,o)i = i^ r ^-i= 



2i x/4^rt 



4f 



St 



where we have used the relation ye v <e y l 2 . Hence (|3.2ip simply reduces to 

\G( s (t), t, e, o) - G(s(t), t, e, o)i < ^o(m, 2*, e, o))k*) - *(t)i. 



Integrating in £, together with (|3.18p yields 

A < C|6|||«j||||M-M|| 



(l + 2t) 3 / 2 - 1 
iV2 



Since limt" 1/2 ((l + 2i) 3/2 - 1) = 0, for a sufficiently small we have Ai < \\\M - M\\. To 
estimate we proceed as follows 



|A| <2 



M( T )G x (s(t), t, s(t), t) - M{ T )G x (s{t), t, s(t), r)di 



M(t)G x {s(t),t, s(t), t) - M(t)G x (a(t), t, s(r), r)di 



= : Ai + A 2 - 
The Lipschitz bound (|3.11j) for s yields 

\n( a (*\ + .u\ V<- 1 l s (*)~ s ( T )l <r (N + 2m 
|Ga;(s(i),i,s(r),r)| < — =— ttttt- < 



2\/4^ (t-r) 3 / 2 



2^ (t-r)V2' 



and consequently 



|Ai| <2||M-M|| / \G x (s(t),t,s(T),r)\dT 



<C\\M-M\\ 



t (t_ 1 T )i/2 dT - Cm||M _ a7||(t1/2 - i^ l|M _ ^ 
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for <7 small enough. To estimate A22 proceed as follows: 
\G x (s(t),t, s(r), r) - G x (s(t), *> 5(r), t)| 
= C 



S(<) S(T) G( S (t), t, S (r), r) - S{t l HT) G(s(tU S(t),t) 



< C 



t - T 

s(t) - s(t) s(t) - s(t) 



+ C 



t - T 

s(t) - S(t) 



t - T 

G(s(t),t,s(r),T) 



t - T 



t - T 

\G(s(t),t,s(T),T)-G(s(t),t,s(r),T)\ 



In order to estimate B\ we use the mean value theorem 

[s(t)-s(t)}-[s(r)-s(T)} 



t - T 



= s(r) - s(t) 



(3.22) 



for some < r < t. By the previous equality and (|3.19[) we have 



Bi < 



C 



:\s(f) ~ s(f)| < 



C 



\M-M\\. 



( t _ T y/2\°vj °VJ\ ^ ( t _ r )l/2 

On the other hand, to handle the term B2, we first note that 

\G{s(t),t,s(T),T) - G(s(t),t, s(t),t)\ 

-(^)-.?(r)) 2 + ( S (t)- S (r)) 2 



<G(s(t),t, S (T),T) 



1 — exp 



(3.23) 



4(t - t) 
Define now 

S := { S {t)-s{r)f -{~s(t)-~ S {r)) 2 = [*(*) - s(r) + g(t) - S(r)} [a(t)-S(t)-{s(r)-S(r))]. (3.24) 
The mean value theorem (|3.22l) and the estimate (|3.19[) lead to 

\[s(t)-~s(t)]-Hr)-~s(r)]\ = |i(f)-T»(f)|(t-r)<C||M-M||(t-r). 
On the other hand, we recall again the Lipschitz estimate (|3.20[) . i.e., 

max{|s(i) - s(t)\, \s(t) - s(r)|} < Cm(t - r), 
for a constant depending on |6o|, which yields an estimate for p.24p . 

\S\ 



(3.25) 
(3.26) 



t - T 



< Cma\\M-M\\ 



The combination of the above inequality with f|3 . 23[) together with the mean value theorem 
shows that 

\G(s(t), t, s(t), t) - G(s(t)A ~s(t),t)\ < G(s(t),t, s(r), T)Cma\\M - M\\, 
and thus the term B2 is estimated using (|3.26l) 

B 2 <Cm^M-M\\a-^-^. 

Multiplying B\ + B2 by M(t) and integrating over the interval [0,i] yields 

A 2 2<Cm f \G x (s(t),t,s(r),T) - G x (s(t),t,s(r),T)\dT 
Jo 

rt 1 
<Cm (B 1 +B 2 )dT<Cm 3 \\M-M\\cr 1 / 2 < —\\M - M\\, 

.In 12 
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for a small enough. 

The next step is to estimate Az- Split the integral into two terms 



M(T)G x (s(t), t, ai (r),r) - M(r)G x (s(t), t, Sl (r), r)dr 



M( T )G x {s(t),t, «i(t),t) - M{T)G x (s{t), t, si(r), r)d 7 



-^31 + -^32- 



The estimate for A31 is very similar to that of J3 from (|3.15p . Indeed, we can analogously 
deduce t 

|^3l|<2||M-M|| / \G x { S (t),t, Sl (T),T)\dT 

Jo 

/°° 1 2 1 

-e- 7 - 2 dz < —\\M-M\\, 
a z 12 

■TSa 

where we used that s(t) — Si(t) > A > for <r sufficiently small with 

A fM-|&ok for6<0 
\M - (M +m)ff for6>0' 

To bound A32 we split 

\G x (s(t),t, Sl (r), r) - G x (5(t), f, Si(r), r)| 
«(t)-«i(r) 



= C 
< C 



G( S (t),t lSl (r),r) - S(t \ SliT) G(s(t),t,~si(r),r) 



t - T 

s(t)- Sl (r) s(t)-si(r) 



t - r 
5(t)-5i(r) 



t - r 



t - T 

G(s(t), t,fli(r),r) 



t-r 
B[+B' 2 . 



\G(s(t), t, si (t), t) - G(s(t), t, h(r), T )\ 



We observe that B[ is estimated exactly the same way as B\. This is a consequence of 

[s(t) - 5(f)] ~ [si(r) - Si(r)] = [*(*) - 5(t)] - [s(r) - g(r)]. 
We can continue from (|3 . 22[) as before to obtain 

C 



\B[\< 



(t-r) 1 / 2 



\M — M\\. 



The estimate for B' 2 is slightly more involved. We write 

B'o = C m ~ gl(r) \G(s(t),t,h(r),T)\ fl-exp f 



t - T 



4(t - r) 



(3.27) 



(3.28) 



for 



S': = -( S (t)- Sl (r)r + (~s(t)-h(r)) 2 

= [*(*) - h(r) + s(t) - Sl (r)] [S(t) - h(r) - s(t) + Si (t)]. 
By the definitions of Si and Si (see (|2.2p ) we have that 

\s(t) - h(T) - s{t) + Si {t)\ = \s(t)-s(T)-s(t) + s(r)\<C\\M-M\\{t-T), (3.29) 

where in the last inequality we have used estimate (|3.25|) . On the other hand. 

\[s(t) - si (t) + s(t) - Si (t)}\ < \s(t) - s(r)| + \s(t) - s(t)| + 2|vjt|V2f + T < Cma, (3.30) 
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if we use again the Lipschitz estimate (|3 . 26[) . Hence putting together (|3.29[) and (|3.30[) we get 
again that 

I 9' I 

-!— L < Cma\\M-M\\, 

t-T - " " 

and consequently, (|3.28[) reduces to 

\B> 2 \ < Cm S ^' ~ Sl{j > G{s(t),t, S!(t), r)a\\M - M\\. 

t — T 

Integrating the previous expression, using the inequality ye~ v < e~ v I 2 , and noting that s(t) — 
Si(r) > A > 0, we can give a very rough estimate that is enough to our purposes: 

/ \B' 2 \dT < Cma\\M- M\\. (3.31) 
Jo 

Thus, from the estimates for B[ and B' 2 from (|3.27[) and (|3.31[) respectively, 

\A 32 \ <Cm [ (B[ + B' 2 ) dr < C\\M - M\\ (ma 1 ' 2 + m 2 a) < i||M - M||, 
Jo V / 12 

for some suitable a small enough. Then, adding the estimates obtained for A\ , A2 and A3 yields 
that T is a contraction satisfying for some a small enough inversely proportional to m: 

\\TM-TM\\ < -\\M-M\\. 
This concludes the proof of Theorem 13.21 as desired. □ 



3.3 Recovery of u 

Theorem 13.21 shows that we have short time existence of a mild solution for problem (|2.2j) (i.e., 
a solution in the integral sense). However, one can easily show that: 

Corollary 3.3. There exists a unique solution of problem (|2.2|) in the sense of Definition [2~2l 
forte [0,T]. 

Proof. Once M is known, one can construct u from Duhamel's formula p. 61) . The smoothness 
and decay of u follow immediately from here. One needs to check also that u has well defined 
side derivatives at s\. But this follows from formula (|3.7|) applied to s\(t) and the estimate for 
\G x (s\(t),t, si(t),t)\ that follows similarly as the calculation in (|3.14j) . □ 

This completes the proof of Theorem 13. II 



4 Proofs of the Main Results 

From the previous arguments, see p.lOj) . it is clear that the obstacle for long time existence in 
this case is the possible blow up in time of t)\\ particularly at the free boundary, i.e. the 

blow up of M(t). We now formalize this idea by showing that we can extend the solution as 
long as the firing rate M(t) is bounded. 

Proposition 4.1. Suppose that the hypotheses of Theorem 13.11 hold and that (u(t), s(t)) is a 
solution to (12.21) in the time interval [0,T]. Assume, in addition, that 

Uq := sup \u x (x,to — s)\ < 00 and that M* = sup M(t) < 00 , 

i£ ( — oo,s(t — e)] *G(*o — £,*o) 

for some < e < to < T. Then 

sup {Im^cc, t)\ with x s (—00, s(t)] , t G [to — e, to)} < 00 , 
with a bound depending only on the quantities M* and Uq- 
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Proof. Differentiating (|3.6[) in x yields 



/s(to-e) 
G(x, t, £, t - e)ua:(£) *o - e)d£ 
-00 

M(T)Gf a ,(s ) t,«(r),r)dr+ / M{r)G x {x,t, ai(r),r)dr 
=:/i-/ 2 +l3. 

The estimate for ^ is straightforward from heat kernel properties and it depends only on Uq- 
Let us deal with I 2 ■ Since M is uniformly bounded in the whole interval to — e < t < t , we get 

\h\<C f \G x (x,t,s(r),T)\dT. (4.1) 



Next, it is shown in [7J Eq. (1.16), pag. 219] that for any Lipschitz continuous function s, 
/■* \ x — s(t)\ 

/ ' v ' l G(x,t,s(T),T)dT < C, te(*o-e.*o), 

for some C depending on the Lipschitz constant of s, to and e. Then this formula allows to 
control the expression in (|4. 1 1) in order to bound the term Ii for t G (to — £, to)- However, this 
bound may depend on to and M* since the Lipschitz constant of s does, see (|3.17|) . 

Finally, the same argument works for ^3, replacing s by s\ in the previous calculations. □ 

With this result in hand, our solutions can be extended to a maximal time of existence and, 
we can characterize this maximal time. The following result holds no matter the sign of b. 

Theorem 4.2. Suppose that the hypotheses of Theorem 13.11 hold. Then the solution u can be 
extended up to a maximal time < T* < 00 given by 

T* = sup{t > : M(t) < 00} . 

Proof. Assume that the maximal time of existence of a classical solution (it(t), s(t)) to (|2.2|) in 
the sense of Definition 12.21 is T* < 00, if not there is nothing to show. By definition we have 
T* < sup{t > : M(t) < 00}. Let us show the equality by contradiction. Let us assume that 
T* < sup{t > : M(t) < 00} and then, there exists < e < T* such that 

M* = sup M (t) < 00 . 

te(T*-e,T*) 

Let Uo be defined as in Proposition 14. 1 1 with to = T*. Applying Proposition ^. 1[ we deduce that 
u x (x, t) is also uniformly bounded for x <G (—00, s(t)] and t e [T* —e, T*) by a constant, denoted 
U* . The same Proposition tells us that U* only depends on M* and also on Uo, i.e., the uniform 
bound of u x (x,T* — e) for x G (— oo,s(T* — e)]. Therefore, we can now restart by using the 
local in time existence Theorem 13 . 1 1 using as initial time any to G [T* — e, T*) for a time interval 
whose length does only depend on U*. Thus, we can extend the solution (u(t),s(t)) to (|2.2[) 
after T* and find a continuous extension of M(t) past T*. Wc have reached a contradiction, 
hence the conclusion of the Theorem follows. □ 

We now show, following Friedman's ideas [7], that it is possible to extend the solution for a 
short (but uniform) time e for b < 0. 

Proposition 4.3. Suppose that the hypotheses of Theorem 13.11 hold and that (u(t),s(t)) is a 
solution to (|2.2|) in the time interval [0,io) for b < 0. There exists e > small enough such 
that, if 

sup \u x (x,t - e)\ < 00, (4.2) 

x£( — oo,s(£ — s)] 
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for < e < to then 

sup M(t) < oo. 

to—e<t<to 

Although the estimate depends on the bound (|4.2[) . e does not depend on to- 



Proof. We use the integral formulation (|3.8I) for M, this time with initial condition at time to — £ 
for some fixed e chosen below, and t E (to — £, to). Thus 

M(t)=-2 { G(a(t),t,£,to -£)«*(£, *o-e)d£ 

*/ — OO 

+ 2/ M(T)G x ( S (t),t,s(T),T)dT-2 [ M(T)G x (s(t),t, Sl (T),T)dT 
•/to— e Jt — £ 

= : K x + K 2 + K 3 . 

Since s(t) > s(t), it follows that G x (s(t),t, s(r),r) < 0. Moreover, M > 0, hence < and 
this term can be discarded. To estimate K3 let 

<E>(i) := sup M(t). 

t -£<T<i 

Note that 

\K 3 \<*(t) [ \G x (s(t),t, Sl (T),T)\dT. (4.4) 

To estimate the derivative |Ga;(s(t), t, s±(t), t)\ we use that the nonlinear part of s is an increas- 
ing function in the case b < as in (|3.13p . Thus, for e small enough, we conclude that 

s(t) - «i (t) = «(t) - s(r) - v r cc- x (t) > \v R \ - \bo\e > 0. (4.5) 

Hence, we can recall the computations in p.!4[) to estimate 

/ |G a (8(t),t, fl i(r),r)|dr<C/ -e"^ dz < C -e~ z " dz. 



It is clear that this last integral can be made less than 1/2 for some small enough e, which is 
independent of to. Substituting the above inequality into f|4.4[) we have the estimate \K$\ < 
i<f>(t). Finally, note that \K±\ < C depending on sup \u x (x, to — e)\ Combining the estimates for 
K t , K 2 ,K 3 with (|4T3j) yields 

M(t) < h(t) + c. 

Taking the supremum on the left hand side, we get that <£(£) < 2C, for all t E (to — £,to), as 
desired. □ 

Remark. Let us point out that the key estimate (|4.5|) to get the uniformity of the time interval 
with respect to to comes from the fact that the nonlinear part of the free boundary s(t) is 
monotone increasing in the case b < 0. For the case b > 0, instead of (|4.5[) we got (|3.16[) . which 
makes impossible to get a uniform estimate with respect to to since m will depend on it. 



Combining Proposition ^. 3[ Thcorcm l4.2[ and 14. II gives global existence for b < 0, as summarized 
in the following result: 

Theorem 4.4. Let uj(x) be a non-negative C 1 ((— 00, s/]) function such that ui(si) = 0, and 
Ui,(ui) x decay at —00. Then there exists a unique global classical solution (u(x,t), s(t)) of the 
equation (|2.2p with b < in the sense of Definition al?^ with initial data uj. Furthermore, the 
function s(t) is a monotone increasing function of t if both b and bo are negative. 

With this the proof of our main Theorem 1 1.1 1 is complete. We emphasize that our Theorem 
321 characterizes the possible blow-up of classical solutions in finite time as the time of divergence 
of the firing rate N(t). 
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5 Study of the spectrum 

In this section we study the spectrum of the linear version /j, = of (|1.6p : 

p t - d v (vp) - d vv p — N(t)8 v = VR on (-00, 0), 

where 

N(t) = -p v (0,t), p(0,t) = 0. 
The objective is to solve the eigenvalue problem 

J d vv p + d v (vp) - Pv(0)5 v = VR = Xp, v G (-00, 0), 
\p(0) = 0, 

with eigenfunctions in the space Ll xp (M.) defined as 

L 2 exp (R) := [p G L 2 (R) : ||p|| iLp(R) < 00} , 

with norm 

2 



(5.1) 



Note that although problem (|5.1[) is only defined in (—00, 0), it can be easily extended to R by 
odd reflection. Following an idea developed in [111 I12j , we consider the equivalent problem to 
(|5.ip defined as 

d vv p\ + d v {vp x ) = Xpx in (-oo,i> H ) U (ur,0), (5.2) 
with p> satisfying the following properties: 
(Fl) PA G LLpW, 
(F2) p A (0) = 0, 

(F3) Matching condition: P\(Vr) = P\(v^), 
(F4) Jump condition: d v p\{v^) = d v p\(v^) + d v p\(0). 
We are going to define the solution for (|5.ip of the form: 



P\(v) = X(-oo,vr-)P L (v) + X(v R ,o)P 2 (v), (5-3) 

where each p l (v), i = 1,2, is a linear combination of the two linearly independent solutions of 
(|5.2p in R, and such that the combination (|5.3[) satisfies (F1)-(F4). 

The functions and p 2 (v) will be calculated by a standard classical method used to 

compute the spectrum for the classical Fokker-Planck equation given by C(p) = Xp, v G R, with 

C{p) := d vv p + d v (vp). (5.4) 

Define first the eigenspace 



C 2 n (R) = { P eL 2 (R) : \\p\\ Li 



with norm 



= / (H-v 2 ) m |p(t,)| 2 di;. 



For completeness we recall a well known result on the spectrum for the classical operator £, 
see for instance [91 [2T] : 
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Lemma 5.1. For any m > 0, the spectrum of the operator C defined in (|5 .4|) on Lf n (M.) is given 
by 

a(C) = {A e C : 9t(A) < § - m} U {-n : n e N U {0}} . 

Moreover, if m > ^ and if n G N U {0} satisfies n + ^ < m, then X n = —n is an isolated 
eigenvalue of C, with multiplicity one, and eigenfunction given by the n-th Hermite polynomial 

dV 
dv r > 



H n (v) = (-1)V» 2 / 2 — e 



v 2 /2 



exp ( 



reduces to the 



In particular, the spectrum of the Fokker-Planck operator L in the space L 
eigenvalues 

\ = -n, ?i€NU{0}. 

It will be very illustrative to give a sketch of the proof of the above Lemma in view of the 
computations for (|5.1D . 



Proof of Lemma \5.1[ Given m <G N, we seek a solution p £ L? m (M) for 

Pw + dv{vp) — Xp, Ae C. 
The Fourier transform of a function p is defined as 

1 



(5.5) 



p(0 



'2it 



p{v)e 



dv. 



Then, the Fourier transform of (|5.5[) yields the following first order differential equation 
which has solutions given by 



m 



for £ > 0, 



(5.6) 



for some constants a, ft £ C. By inverse Fourier transform we get: 
P(«) = 4= / e^p(0dC 

V Z7T Jm 



«£-A e -? 2 /2 



d^ + p 



(5.7) 



3 w£ 



(-0 



This is, for each A G C, p is the linear combination of the two linearly independent solutions of 
(|5.5[) . The constants a,/3,\ will be determined from the boundary conditions. 

Note that the Fourier transform is an isomorphism from L^(R) to H m (M). Thus, A belongs 
to the spectrum of C if and only if the function p from (|5.6[) belongs to H m (M.). Since p is 
sufficiently smooth and rapidly decaying for any m, provided we stay away from the origin, we 
just need to check if p £ H m (0), i.e., all derivatives of p of order less or equal than m arc square 
integrable near the origin. 

Clearly, the values A = — n, neNU {0} are special and will be considered later. In order to 
see if any other values of A G C are admissible, we compute a general n-derivative of the term 



£"/ 2 around £ = : 



It is easy to see that 



Q{n) 



-> p -« 2 /2 



£ L 



-A— n 
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if and only if Re(X) < — m. Since the spectrum of £ is closed, this shows that 

a{C) D {A G C 31(A) < \ - m}, 

as claimed. 

Now we study the values A = —n, n = 1, . . . , m — 1. Let p„ be the corresponding eigen- 
function, with Fourier transform given by (|5.6[) . Since p n belongs to H m near the origin, its 
m-derivative will be L 2 integrable, while the rest of the derivatives p n , d^p n , . . . , dj?" (p n ) 
will be continuous at the origin. This forces to have a very precise values for a, ft. In partic- 
ular, for n even, a = /?, while for n odd we must have a = —(3. Then we have shown that 
A = —n, n = 1, . . . , m — 1 are admissible eigenvalues and the corresponding cigenspaces are one 
dimensional, with eigcnfunctions given by the well known Hermite polynomial H n as 

p n (v) = e- v2 ' 2 H n {v). 

□ 



We consider now the original problem (|5.1[) and seek for solutions p(v) of the form (|5 .3|) . 
Our hrst observation is that the values for A are determined only by the decay of p as v — > — oo. 
Consequently, if we impose that the function p 1 belongs to L 2 xp (M.), then this fixes the possible 
values of the eigenvalues A as in Lemma 15.11 In particular, there is no continuous spectrum. 
Moreover, for each A„ = —n, n £ N, we must have 

p\v) = aH n {v)e- v2 '\ (5.8) 

for some ael. 

The difference between our problem and the classical Fokker-Planck operator lies in the fact 
in the interval (vr, 0) all solutions to the ODE (|5.5[) for A = — n as given in (|5.7p are admissible 
since the behavior at infinity does no play any role. Nevertheless, we can find a better expression 
for the two linearly independent solutions in ()5.7[) . One of those is the well known (|5.8|) . The 
other solution can be easily found by making the following ansatz: 

p\v) = e~ v2 / 2 H n (v)g(v). 

By imposing that p 2 (v) satisfies (|5.2|) one can obtain an equation for g(v) that reads as 
follows: 

2H'J - vg'H n + H n g" = 0. 
This equation has the following general solution: 

f v e s2/2 

9{v)= ^Lws) ds+ ^ 

for some constants /3i, /?2 S K, and where we have fixed any vq £ (vr, 0) such that H„ (vq) ^ for 
the integral to be well defined. Note that g is well defined for all v even where the denominator 
vanishes because the Hermite polynomials only have single roots. Consequently we define 

P 2 (v) :=f3 ie - v2 / 2 H n (v) ^ -^-ds + /3 2 e- v2 / 2 H n (v), 

and the eigenfunction corresponding to A = — n is simply 

ae~ v 1 H n (v), i)6 (—oo,vr), 

Pn(v) = ( 2 , r p^ 2 /2 2; (5.9) 

f3 ie ^ 2 / 2 H n (v) / ^^ds + foe- v2 / 2 H n (v), v £ (v R ,0]. 
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for some real constants a, (3 For simplicity, let 



e ,'/2 



It is clear, by doing a careful Taylor expansion, that if v± is a root of H n , then there exists a 
finite limit for 

A Vl ,n '■= lim 6 n (v) ^ 0. 
v— yv\ 

Now we are ready to check if (|5.9[) is an admissible eigenfunction. In the case n is odd 
integer, the Hermite polynomial H 2n+ i vanishes at zero, but as we have mentioned, 

02n+i(v) -> Ao, 2ll +i 7^ 0, as u -> for any n e N. 

Then in this case condition (F2) is satisfied only when (3i = 0. Then, if we wish p 2n +i to be 
a continuous function as stated in condition (F3), we must have a = (3 2 unless H 2n +i(vR) = 
that will be considered afterwards. The solution constructed this way does not satisfy condition 
(F4) , so we conclude that 2n + 1 is not an admissible eigenvalue. 

On the other hand, let us check if p 2n is an admissible eigenvalue. For even integers it holds 
that H 2n (0) ^ 0. Thus we can simply take vo = 0. Consequently condition (F2) is satisfied if 
and only if 

13 2 = 0. (5.10) 

The matching condition (F3) implies 



aH 2n (v R ) = P\H2 n {v R ) I ds. (5-11) 

JO H n( s ) 

Here we do need to distinguish two cases: if Vr is not a root for any H 2n , then the above equality 
implies 



Q=A /. Tm d '- (5 - 12) 

If instead H 2n (vn) = (note that Hermite polynomials only have single roots), one can repeat 
a Taylor expansion around vr for 9 2n (v) and see that 

2n {v) -> A VRt2n ^ 0, as v -> vr. 

Consequently (|5.11[) cannot be satisfied for these n such that H 2ti {vr) = 0. 
Using conditions (|5.10|) and (|5.12[) for p n we get 

r v R e s 2 /2 

— 2-prds, V S (-QO,Vr), 

„ e 2 ;; 2 (5.i3) 

— 2—-ds, se(«jj,0]. 

'0 U 2n\ S ) 

One can easily check that the jump condition (F4) is satisfied if and only if 

H 2n (0) =H 2n (v R ). (5.14) 

Summarizing, we have shown the following: 
Theorem 5.2. Consider the operator 

j d vv p + d v (vp) - p v (0)A v=VR = Xp, v 6 (-oo, 0] 

\p(o) = o, 

subject to conditions (Fl) - (F4). 
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1. There is no continuous spectrum. 

2. The value A = is an eigenvalue with a one- dimensional eigenspace spanned by the function 



Poo{v) 

for 



e " 2/2 v E (-oo,v R ), 

a e~ v2/2 / e v2 ' 2 dv v e (v R ,0], 



I e" 2/2 dv 



I VR 

3. There exists a countable set S C R such that for all vr S, there are no other eigenvalues. 

4- If n and Vr happen to satisfy the compatibility condition (|5.14[) , then \ = —In is an 
eigenvalue with eigenspace of finite dimension spanned by the eigenfunction pi n (v) defined 
in flUBJ) . 

Remark. We remark that the steady state p^ (v) was previously obtained in [TTJ [3] . In this 
last paper, it was also shown exponential decay towards equilibrium p^. However, the speed of 
convergence is unknown and the spectral analysis does not seem to give any insight. 
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